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Abstract 



^ \ The entrainment matrix (also termed the Andreev-Bashkin matrix or the mass- 

' density matrix) for a neutron-proton mixture is derived at a finite temperature in a 

neutron star core. The calculation is performed in the frame of the Landau Fermi- 

OO 

\ the temperature dependence of the entrainment matrix is described by a universal 

■ function independent on an actual model of nucleon-nucleon interaction employed. 

. The results are presented in the form convenient for their practical use. The en- 

trainment matrix is important, e.g., in kinetics of superfluid nucleon mixtures or in 
O ' studies of the dynamical evolution of neutron stars (in particular, in the studies of 

star pulsations and pulsar glitches). 
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1 Introduction 



It is well known, that the neutron star core becomes superfluid (supercon- 
ducting) at a certain stage of neutron star thermal evolution (see, e.g., Ref. 
[1]). It is generally agreed that protons pair in the spin singlet {^Sq) state, 
while neutrons pair in the spin triplet (^^2) state in the neutron star core. A 
variety of different models of nucleon pairing have been proposed in literature 
(references to original papers can be found in Yakovlev et al. [2] and in Lom- 
bardo and Schulze [1]). These models predict very different density profiles 
of neutron (n) or proton (p) critical temperatures Tcn,p(p). In addition, it is 
not absolutely clear whether the actual projection of angular momentum of 
neutron pair onto quantization axis is rrij = (as one usually assumes) or it 
can be \mj\ = 1 or 2. For example, Amundsen and 0stgaard [3] found that 
the energetically preferable state of a neutron pair can be a superposition of 
states with different mj. 

Despite many theoretical uncertainties it is obvious that superfiuidity strongly 
affects the evolution of neutron stars, for example, its coohng (see, e.g., Ref. 
[4]), neutron star pulsations (see, e.g., Refs. [5] - [10]) and is probably related 
to pulsar glitches (see Refs. [11,12]). 

One of the key ingredients of hydrodynamics and kinetics of superfiuid mix- 
tures is the entminment matrix p^a' (also termed the Andreev-Bashkin matrix 
or the mass-density matrix). Implying, for simplicity, that the only baryons in 
a neutron star core are neutrons and protons, the matrix p^a' can be defined 
as [13]: 



"n \Pn Pmi qp ~r Pnn ns + Pnp Vps , (1) 

•^P = (Pp - PPP - Ppn) ^qp + Ppp ^ps + Ppn V^s ■ (2) 

Here pa — TnaUa] ria and nia are the number density and the mass of nucleon 
species a = n or p; and Vas are the mass current density and the superfluid 
velocity; Vqp is the normal velocity of thermal excitations (see, e.g., Refs. 
[14,15]). We assume that Vqp is the same for nucleons of both species. Eqs. 
(1) and (2) differ from a "natural" expression for the mass current density 
Ja — PaVa (with Vq being the momentum per unit mass of nucleon species 
a) for two reasons. 

First, three independent motions can exist in a mixture of two superfiuids, 
each carrying a mass (see, e.g., Ref. [16]). They are the motion of thermal ex- 
citations with the velocity Vqp and two superfiuid motions with the velocities 
Vns and Vpg. 

Second, the superfiuid fiow of one component of the mixture entrains a fiow 
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of another component and vice versa. For example, the superfluid motion of 
neutrons carries along some part of the mass of protons because nucleon liquids 
are strongly interacting (see, e.g., Ref. [17]). The non-diagonal elements pnp 
and Ppn, therefore, characterize the intensity of neutron-proton coupling. In 
particular, in the absence of interactions between neutrons and protons one 

has Pnp = Ppn = 0. 

It follows from the phenomenological analysis of Andreev and Bashkin [13] 
that the entrainment matrix should be symmetric (pnp = Ppn) , since it can be 
presented in the form: 



where E is the energy density in the coordinate frame in which Vqp = 0, and 
T is the temperature. 

The entrainment matrix p^a' of a neutron-proton mixture was calculated by 
Borumand et al. [18] and Comer and Joynt [19] for T — 0. However, in many 
cases the zero-temperature approximation cannot be justified. For example, 
one needs the matrix paa' at non-zero temperatures for analyzing kinetic prop- 
erties of matter, especially the kinetic coefficients (the bulk and shear viscosity, 
the diffusion coefficient). Also, the entrainment matrix is required at T for 
investigating pulsations of warm neutron stars or stars possessing the pulsa- 
tion energy of the order of or higher than its thermal energy (It is implied that 
pulsation energy, being dissipated, is able to heat the star substantially and 
thus to change Paa'- A simple example, illustrating the influence of tempera- 
ture related effects on the neutron star pulsations, was considered by Gusakov 
et al. [20]). 

In this paper the entrainment matrix of a neutron-proton mixture is derived for 
non-zero temperatures. Calculations are performed in the frame of the Landau 
Fermi-liquid theory, generahzed by Larkin and Migdal [21] and Leggett [22] to 
the case of superfluidity. For the sake of simphcity we assume the singlet-state 
(^5*0) pairing of nucleons of both species. In Section 4.1 we will show how 
the obtained results can be extended to the case of triplet-state neutron 
pairing. 



2 A current-free neutron-proton mixture 

Before calculating the matrix Paa', let us consider a neutron-proton mixture 
in the absence of currents. A simple generalization of the Hamiltonian of a 
superfluid Fermi-liquid, suggested by Leggett [22], to the case of superfluid 




(3) 
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mixtures, gives 



H - /i„Nn - //pNp = Hlf + Hpairing ■ (4) 

Here, H is the hamiltonian of the system which is the sum of the stan- 
dard Fermi-hquid hamiltonian Hlf for mixtures and the pairing hamiltonian 
Hpairingi and Ha are the number density operator and the chemical poten- 
tial of nucleon species a, respectively. The expression for Hlf has the form 
(see, e.g., Ref. [17]): 



pp'aa'aa' 

In Eq. (5) the summation is taken over the particle momenta p andp', as well 
as over the spin projections cr and a' onto the quantization axis and over the 
particle species a,a' = n or p; = = or is the annihilation 
operator of a quasiparticle (not the Bogoliubov excitation!) of species a in a 
state (per). We restrict ourselves to a spin-unpolarized nucleon matter. This 
allows us to simplify the notations. For instance, we will drop spin indices, 
whenever possible. Furthermore, ~ ^ (Ppa ~ li'D' where 9{x) is the step 
function; Sq"^ (p) = "^Fadpl — Pfq)' where v-p^ and are, respectively, the 
Fermi- velocity and Fermi-momentum; {p,p') is the spin-averaged Landau 
quasiparticle interaction (we disregard the spin-dependence of this interaction 
since it does not affect our results); /"°'(p,p') is the (spin-averaged) second 
variational derivative of the energy with respect to the number of particles. 
Therefore, it is invariant under transformations p ^p' and a ^ a' (see, e.g., 
Refs. [14,17]): 



n' ip,p') = (P,P') = n' (P',P) ■ (6) 
The pairing hamiltonian can be written as: 

HpaiHng = E V^"^ (P,P') a^l' aL;l\ aL;\ a^^ . (7) 

pp'a 

Here we assume the following symmetry conditions: 

V(°) (p,p') = V(") (p',p) = V(°) (p, -p') = V(") (-p,p') . (8) 

If the matrix element V*-"^ (p,p') does not satisfy these conditions (e.g., for 
singlet-state pairing of nucleon pairs with non-zero orbital angular momen- 
tum), it should be symmetrized in such a way as to obey Eq. (8) (see, e.g.. 
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Ref. [23]). To diagonalizc the operator (4) we introduce the standard operators 
of Bogohubov excitations : 



«ff = + (9) 

= <^^£^-^^^^^-;r (10) 

The parameters u^"^ and t;^"^ are related by the inversion symmetry and nor- 
mahzation 

<^=«i\ v^^^=v^^i, = (11) 

Now one can find the eigenvalues of the operator (4) 



prra 

+ ^V(")(p,y) (l - 2f(r)) (l - 2f^)) , (12) 

pp'a 

where and f^") are the distribution functions of quasiparticles and Bogoh- 
ubov excitations, respectively, given by 



n^^^ = Hfalfl) = (|aif aif I) = + - t\ (13) 

r= (l^if ^ifl) = (l^'i^'^'i^D- (14) 



The entropy of the system is given by the usual combinatorial expression: 



(15) 



p<Ta 



Minimizing the thermodynamical potential F = E — finnan — fJ^pUp — TS with 
respect to unknown functions f^") and {v^"^ can be expressed through 
using Eq. (11)), one finds 



^/I(-^Hp}T^, (16) 

(17) 

Here and below we use the system of units in which h = = V = 1, where h is 
the Planck constant, k-Q is the Boltzmann constant, and V is the normalization 
volume. In Eqs. (16) - (17) is the energy of a Bogoliubov excitation with 



u. 



(«)2 
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momentum p. The superfluid gap A^'*) of nucleon species a can be determined 
from the equation 

A(") = -Y^ v(") (p,p') u^fv^^ (l - 2fJ^)) . (18) 
p' 

Finally, e*^") (p) is the quantity which formally coincides with the energy of 
quasiparticle species a in the mixture of non-superfiuid Fermi-liquids, 

(p) = (p) + E M) {n^'^ - ^''0 ■ (19) 

p'a'a' 

In Eq. (19) the quasiparticle distribution function n'^^ is determined by Eq. 
(13) with f^°) and u'^^ taken from Eqs. (16)-(17). The first term in the right- 
hand side of Eq. (19) can be estimated (p) ~ (t + A(")), where A^") 
is a typical value of the gap. In thermodynamic equilibrium, the second term 
in Eq. (19) is much smaller than the first one because for any function /(p), 
smooth in the vicinity of a Fermi surface, one has the estimate: 

r m P' {n(^^ - Oi"'^) dp ~ /(PfJ {[T/^iX + [A^'^V/^a]^) , (20) 

J 

Thus, since <S 1, the second term in Eq. (19) can be neglected. 



3 A neutron-proton mixture with superfluid currents 

3.1 General consideration 

In a system with superfluid currents the plane-wave states of nucleons (p + 
Qq,, T) ^-iid (— p+<9a' i) paired (note, that we consider singlet-state pairing 
of both species and assume = rUaVas <^ Ppa)- ^^^^ pairing 
hamiltonian should be written as (see, e.g., Refs. [24,25]) 

Hpairing(<3a) = E ^^P') '^P'+QJ^^-pUq^I^-P+Q^I^^^ (21) 

pp'a 

Here Vq*"* (p,p') is the matrix element for the scattering of a pair of quasipar- 
ticles species a from states (p + Q^, ]), (— p + Q^, J,) to states (p' + Q^, T); 
(— p' -|- Qq,, i). Due to the rotational invariance, the expansion of Vq"'' (p,p') 
in powers of will contain the terms ~ and higher. Therefore, as we 
will work in the linear approximation in Q^, we will neglect the dependence 
of v'q' (p,p') on momentum and put Vq^ (p,p') ~ V^°^ (p,p')- 
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Expressing the quasiparticle operators in terms of Bogoliubov excitation op- 
erators 



where t/^"'' and V^""^ satisfy the equahties 

C/^°) = C/5, yp(°) = yi^\ C/^")2 + y;")2 = l, (24) 
similar to (11), one obtains the expression for the energy density: 

pp'tTir'aa' 

pp'o 

X (1 - ^Sq. - ^55h^j (1 - ^u. - ■ (25) 

Here J^^q and J^^pfg are the distribution functions of quasiparticles and 
Bogohubov excitations with momentum (p + Qa), respectively: 

= ^ + ^ ^^^^ - ^;") ^ , (26) 

•^Sq„ = (I«Q.T^'Sq„tI) - i\b';^Qj;-lQj)- (27) 

The entropy of the system is still given by Eq. (15) with the distribution 
function f^") replaced by -^p+Q • The minimization of the thermo dynamical 
potential F — E — ji-aU^ — l-ipUp — TS with respect to q and U^'^'^ yields 



(28) 



-r-(°) _ 
•^p+Qc ~ 


1 


1 + e 


_ 


-'■ / t7-(") t7-(") 

2 l^P+Qa ~ ^-P+Qa 


/-/(°)2_ 


2V 


^P 



+ ^-P+Q. 



-P+Qa P+Qc . 



(30) 



In Eqs. (28)-(30) g^+^g^ is the energy of a Bogoliubov excitation with momen- 
tum (p + Qa)- The stability of the system imphes (Sp+g > 0. Furthermore, 
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is the superfluid gap which can be found from the equation: 

= _ ^ v(") (p,p') t/^fV;-) (i - ^^Q^ - ^L;VqJ . (31) 

p' 

Finally, -f^p+g is the quantity which formally coincides with the energy of 
a quasiparticle with momentum (p + Q^) in the mixture of non-superfluid 
Fermi-liquids, 

p'tr'a' 

(32) 

The quasiparticle distribution function J^^q is defined by Eq. (26). Since 
Qa ^ Ppa' expand -f^p+Q in terms of Qq.' (a' = n, p) and write 

+ (33) 

In the case of singlet-state nucleon pairing there are only three vectors p, 
and Qp which can form the scalar AH^°'\ Neglecting all terms which are 
quadratic and higher order in Qa/Ppa^ write 

A<)=E7aa'(p) (34) 

a' " 

where 7aa'(p) matrix to be derived in the next section on the Fermi 

surface of particle species a {at p = Pp„). Taking into account Eqs. (8), (29)- 
(31), and (33)-(34) and neglecting the terms ~ Q^, one has 

Pi") = A("), f/^") = = vl^'l (35) 

Now the energy of Bogoliubov excitations , as well as the distribution 

functions of quasiparticles and Bogoliubov excitations can be expanded in 
analogy with Eq. (33) as 



'P+Qa 



(a) 



iW-(") 
■^^P+Qa 



n. 



(a) 



+ 



(36) 
(37) 



3.2 The calculation of matrix ^aa' 



To calculate the matrix ■yaa'iPpa) ^^^^ make use of Eq. (32). Restricting 
ourselves to the terms linear in Q^, we expand all functions in Eq. (32) using 
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the formulas (33) and (37). Then, taking into account Eqs. (19) and (20), and 
neglecting all terms in Eq. (32) which depend on (n^? •* — 0^^), we obtain 

Let us calculate the sum in this equation. The main contribution to the sum 
comes from a narrow region of \p'\ ~ Pp^, since the function in the curly 
brackets is essentially non-zero only close to the Fermi surface of particle 
species a'. Furthermore, in a smooth function (p,p') we replace |p| and 
\p'\ by Ppq, and p^^,, respectively, and expand it into Legendre polynomials 
Pi{cos9): 

r"'(p,p')=E/r>Kcose), (39) 
I 

where 9 is the angle between p and p'. Using isotropy of the gaps A^") (a — 
n, p) and Eq. (34) we obtain 



E n' M) ^ ^K'^ - - ^-^^ ^ > (40) 



p'cr' '-'^p 



dEi?'^ ' 3 pp^ 

{a') 



p'a' 



^' " 3 ml, 



(41) 



Here iVoa = m* pp^/vr^; m* = p^^hYa effective mass of particle species 

a. The function which is given by 



*. = -^E^. (42) 



was calculated numerically and approximated by Gnedin and Yakovlev [26]. 
Their fit of is given in Appendix A. When calculating the authors 
adopted the standard approximation in which the dependence of the gap on 
the absolute value of particle momentum is neglected (consequently, in the 
isotropic case the gap is a function of temperature only), A^") pa A'^")(|p| = 
Pp^) = A*^")(T'). In this paper we also use this approximation. 

Now, writing Eq. (38) for neutrons (a = n) and for protons (a = p), taking 
into account Eqs. (34) and (40)-(41), and equating prefactors at the same 
in its left- and right-hand sides, we arrive at the set of linear equations for the 
matrix jaa'- Its solution has the form: 
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m* S 



1 + 



1 + 



013 



171(3 



3 ^mlmp 

TP act 



1 /?V2\ 

S \PFaJ 



3/2 



i^r^ (1 - */3) , 



a0^ 



(43) 
(44) 

(45) 
(46) 



Here a ^ (3 (thus if, e.g., a — n then /? = p). The effective masses m* are 
related to the parameters through the equation (see Refs. [17,18]): 



m„ 3 



(47) 



4 The entrainment matrix 



4..1 Superfluid current and the matrix paa' in different limiting cases 



Eqs. (33) and (36)-(37) allow one to calculate the entrainment matrix p^a' and 
to express it in terms of ^aa'iPhh)- This can be done, for instance, with the 
aid of Eq. (3) and the energy density given by Eq. (25). It is easier, however, 
to obtain Pa^i by calculating the mass current density of quasiparticles 
(we remind that Vqp = in a chosen coordinate frame). When doing so we 
take account of the fact that the expression for the mass current density of 
non-superfiuid Fermi-liquid can be applied to the superfluid state as well (see, 
e.g., Refs. [22,18]). In our case this means that the expression for has the 
form (see, e.g., Refs. [14,27]): 

Substituting Eqs. (33) and (37) into Eq. (48) and performing a simple inte- 
gration, we obtain the formulas, similar to Eqs. (l)-(2), with the entrainment 
matrix equal to (a, a' = n, p): 

Paa' = Pa laaiPpa) ~ " (^^) 

From Eq. (44) it follows that the entrainment matrix is indeed symmetric in 
accordance with Eq. (3), paa' = Pa' a- At T — Tea (where Tea is the critical 
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temperature of quasiparticle species a) one has ^q, = 1, and consequently, 
Paa' = and Ja = 0. 

Let us analyze the matrix p^a' in different limiting cases. 

I. Let Fi"^ = 0, i.e., the interaction between neutrons and protons is absent. 
Then each component a of the mixture can be treated as an independent 
superfluid Fermi-liquid while the matrix p^a' is diagonal. Eqs. (49) and (43)- 
(46) yield 

Pa{l-<l>a) _ _ 

Paa ~ I paa ^3 ' /^np — Ppn — ^ . I^OUJ 

The diagonal elements paa coincide with the well-known expression for the 
superfluid density of one component Fermi-liquid (see, e.g., Refs. [22,28]). 

II. Let the temperature of the mixture be equal to zero. Then = and the 
entrainment matrix can be rewritten in the form 



^Pa^fl + ^l, (51) 



Pnp = Ppn = -^-T ^n^P f^' ' (52) 

in agreement with the results of Borumand et al. [18]. 

III. Finally, let us suppose that the only one component a of the mixture is 
superfluid. In this case $^ = 1 (/3 7^ a) and we have: 



(1 + FrV3) (1 + F^'/i) (rrVs)' 



/ -1 -Tv \ ^^^OL 

Paa — Pa\^ ~ ^a) T / \ / ^ \z 

(1 + Fr/3) (1 + Ff^/3) - (F^^/3) 

= Pnp = Ppn = . 



,(53) 
(54) 



We have calculated the entrainment matrix of a neutron-proton mixture at 
non-zero temperatures assuming singlet-state pairing of nucleons of both species 
and isotropic gaps. However, in reality neutron pairing occurs in the triplet 
state of nucleon pair with an anisotropic gap. In this case the expansion 
(34) becomes invalid because the momentum p is no longer the only vector 
which characterizes the system in the absence of superfluid currents; the quan- 
tization axis specifies an additional direction. Therefore, Eq. (34) should be 
replaced by: 

A//^-)=5^ , (55) 

a' " 

where Gaa' is a matrix composed of vectors. It can be, in principle, obtained 
from Eq. (32) in a manner similar to the derivation of the matrix 7aa'- As 
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follows from Eq. (48), components of the entrainment matrix p^a' will be 
tensors (rather than scalars as in the isotropic case). These tensors will be 
expressed in terms of the Landau parameters F"" , with / > 1. The quantities 
F^°° at Z > 2 are not known for neutron star matter and are not necessarily 
small. Thus, a strict calculation of the entrainment matrix in the case of 
triplet-state neutron pairing is rather complicated (the analogous problem in 
deriving the supcrfluid density was discussed in details by Leggett [23] in the 
context of the anisotropic phase of helium-3). 

To proceed further we follow Baiko et al. [29] and assume that the neutron star 
matter can be treated as a collection of microscopic domains with arbitrary 
orientations of the quantization axis. Then the entrainment matrix, being 
averaged over the domains, will be "isotropic" (its elements will be scalars). 
Thus, one can use Eq. (49) for the averaged entrainment matrix by introducing 
an effective isotropic gap which we choose according to Baiko et al. [29] as 

A2(r)=min{A(-)(|p|=Ppj}. (56) 

Here A^g (T) is obtained as the minimum of the angle-dependent gap A^°^ on 

the neutron Fermi surface. The use of Eq. (49) with the effective gap A^g (T) 
allows one to obtain qualitatively correct results for the matrix p^a' in the case 
of triplet-state neutron pairing. The fit of AiJ(T) for the case of neutron 
pairing with mj = was obtained by Yakovlev and Levenfish [30] and is given 
in Appendix A. 

4-2 Landau parameters 

In order to find the matrix p^a' it is necessary to know the Landau parameters 
F"" (and hence the quantities , see Eq. (46)) for an asymmetric nuclear 
matter. In general, the quantity is a function of the baryon number den- 
sity fih and the asymmetry parameter S = (nn — np)/n\:,. Due to the charge 
symmetry of strong interactions one can write 



/r(nb,5) = /r(nb,-5), (57) 

/r(nb,5) = /r>b,-5). (58) 

The function /""'(nb,5) can be expanded in powers of 5 < 1. Neglecting all 
terms quadratic and higher order in S, Eqs. (57)-(58) yield 

fr = a{nk) + Sb{n^) + 0{S^), (59) 

/r = a(nb)-56(nb) + 0(52), (60) 

/-P = c(nb) + 0(0- (61) 
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This approximation was proposed by Hacnscl [31]. To calculate the functions 
a(nb), &(^b) and c(nb) wc need to know the dependence (ub) for any two 
values of the asymmetry parameter 6 = 6i and 5 = 62 (while the knowl- 
edge of this dependence for more than two values of S would enable one to 
find the terms non-linear in 5 in the expansions (59)-(61); to our best knowl- 
edge, these data are unavailable in the literature). The Landau parameters for 
asymmetric nuclear matter can be calculated microscopically within a nuclear 
many-body theory, starting from the nucleon-nucleon interaction in vacuum. 
However, nearly all existing calculations are limited to a simpler case of a 
pure neutron matter or symmetric nuclear matter. In the case of pure neutron 
matter calculations of the / = 1 Landau parameters were performed in some 
range of neutron matter density above nuclear density (see Refs. [32]-[36]). 
In the case of symmetric nuclear matter, calculations were usually done at 
normal nuclear density (see Refs. [37]-[39]) with a notable exception of the 
calculation of Jackson et al. [33]. We are aware of only one calculations of 
Landau parameters in asymmetric nuclear matter in beta equilibrium with 
electron gas (simplest model of neutron star matter), by Shen et al. [40], who 
however restricted to / = parameters, while we need / = 1 ones. In view of 
this situation, we decided to use the density dependent / — 1 parameters for 
symmetric nuclear and pure neutron matter calculated by Jackson et al. [33] 
more than two decades ago. 

These authors calculated the Landau parameters for symmetric nuclear matter 
{S — 0) and for pure neutron matter {S = 1) using two model potentials of 
nucleon-nucleon interaction: Bethe- Johnson v6 (BJ v6) and Reid v6. In the 
case of the symmetric nuclear matter the quantities /f can be expressed as 



/r(%. 0) = /r(n.. 0) = pm+f;m _ ^^^^ 



2-/Vosym 

/r>..0)^ ^-'"-°'-^-'"^-°^ - (63) 

■^-'^Osym 

Here A^'osym = -/Von(?^b, 0) = Nop{nb, 0). The plots of Fi and F{ versus the wave 
number /cpsym = (37r^nb/2)-'^/^ are given in Figures 14 and 17 of Rcf. [33] for the 
Reid v6 and BJ v6 interactions, respectively. We have fitted these functions 
by simple analytical formulas which are given in Appendix B. Eqs. (59) - (61) 
yield 

«(^b) = /r(nb, 0) = /r(rib, 0) , c(nb) = f^^ny,, 0) . (64) 
Now, considering the pure neutron matter, we obtain 

The function Fi(nb, 1) = -f^i(fcFpure)) with /cppure = (37r^nb)^^^, is plotted in 
Figures 20 and 22 of Ref. [33] for the Reid v6 and BJ v6 interactions, respec- 
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Fig. 1. The quantities Paa' / Pa (a, a' = n or p) versus temperature T for two inter- 
actions: BJ v6 (solid lines) and Reid v6 (long dashes) at the baryon number density 
nb = 3no. The proton (Tcp = 5 x 10^ K) and the neutron (Ten = 6 x 10^ K) critical 
temperatures are marked by the vertical dashed arrows. 



tively. The fitting formulas for -Fi(fcFpure) ^ire given in Appendix B. From Eq. 
(59) we have 

%b) = /r(rib,l)-a(rib). (66) 

Thus, assuming that /f " is linear in 5, one can calculate the Landau param- 
eters of the asymmetric nuclear matter and find the entrainment matrix. 

Our results are illustrated in Fig. 1. We show the temperature dependence of 
Paa' / Pa («, a' = u or p) for the BJ v6 (solid lines) and Reid v6 (long dashes) 
interactions. The number density of baryons is taken to be nb = 3no, where 
no = 0.16 fm~^ is the number density of nuclear matter at saturation. For the 
chosen nb, the equation of state of Heiselberg and Hjorth- Jensen [41] yields 
5 = 0.837. The nucleon critical temperatures are taken to be Ten = 6 x 10^ K, 
Tcp = 5 X 10^ K. This information is sufficient for computing the entrainment 
matrix at any T. However, it should be stressed, that our approach is not self- 
consistent. Strictly speaking, one needs to calculate the equation of state, the 
Landau parameters, and the nucleon critical temperatures using one model of 
strong interactions. 
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At T > Tcp in Fig. 1 the matter is non-superfluid and all the matrix elements 
Paa' = 0. At Ten < T^p protons become superfluid and hence 7^ 0. Fi- 
nally, at T < Ten both protons and neutrons arc superfluid and all the matrix 
elements are non-zero. With decreasing temperature, the elements p^a' rapidly 
approach their asymptotes (51) and (52). Notice, that Ppp/pp depends essen- 
tially on the model of strong interactions employed (the difference between 
the solid and dashed curves marked by pp is large), for Ppn/Pp and Pnp/Pn the 
dependence is even stronger: even the signs of these quantities differs for the 
BJ v6 and Reid v6 interactions (for the Reid v6 interaction ppn = Pnp < 
because < 0; see Eq. (52)). The quantity Ppp/pp at T — > is equal to 
Ppp/Pp ~ 1-19 for the BJ v6 interaction and to Ppp/pp ~ 0.75 for the Reid v6 
interaction. Thus, the value of ppp/ pp is approximately half of that obtained 
in the estimate of Borumand et al. [18], Ppp/pp ~ 2. 



5 Conclusions 

We have derived the entrainment matrix paa' at non-zero temperatures. The 
calculation is performed in the frame of the Landau Fermi-liquid theory gen- 
erahzed by Larkin and Migdal [21] and Leggett [22] to the case of superfluid 
matter. The expressions for p^a' reproduce two limiting cases studied in the 
literature: the case of zero temperature and the case in which the interaction 
between neutrons and protons is absent. 

The results are presented in the form convenient for their practical use. In 
particular, for calculating the entrainment matrix one needs the Landau pa- 
rameters of symmetric nuclear and pure neutron matter. These parameters 
were taken from the paper by Jackson et al. [33] for two model potentials of 
nucleon-nucleon interaction (BJ v6 and Reid v6) and approximated by simple 
analytical formulas. While results for p„„ are quite similar for both nucleon- 
nucleon potentials, the values of ppp differ by some thirty percent, while much 
smaller non-diagonal matrix elements, p„p and pp„ differ in sign. Clearly, re- 
sults for the non-diagonal entrainment matrix elements are very sensitive to 
the nucleon-nucleon interaction. Moreover, as we used old values of the Lan- 
dau parameters, our results should be updated as soon as new, more realistic 
values of the / — 1 parameters become available, obtained, e.g., via modern 
renormalization group approach (see, e.g., Rcf. [36]) applied to pure neutron 
matter, symmetric nuclear matter, and maybe also to the most astrophysically 
interesting case of asymmetric nuclear matter. 

The generalization of the entrainment matrix to the case of three or more 
interacting baryon species is straightforward. It just implies an increased cor- 
responding number of baryon indices in all formulas and an associated increase 
of the dimension of matrices p and 7. 
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The temperature dependence of the matrix paa' , ftiHy described by the universal 
function $(,(T), is known rather rehably. Unfortunately, this cannot be said 
about the values of the Landau parameters. Clearly, new calculations of the 
Landau parameters of asymmetric nuclear matter would be highly desirable. 

The matrix p^a' at non-zero temperatures is needed to study the kinetics of 
the neutron star matter as well as to investigate the dynamical evolution of 
neutron stars, especially, their pulsations. As we have shown, the entrainment 
matrix varies with temperature. Because the matrix paa' enters hydrodynamic 
equations which determine the neutron star pulsations, the frequencies of (su- 
perfluid) pulsation modes should vary with T and hence (due to a star cooling 
or heating), with time. This gives a potentially powerful method to probe 
very subtle properties of superdense matter by measuring the dependence of 
pulsation frequencies on time. We intend to consider the related problems in 
a separate publication. 
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A The analytical approximation of the function 

Eq. (42) for $q can be written as: 



Here v — A(T)/T; A(T) is the gap which depends on the type of pairing. For 
singlet-state pairing of nucleon species a the function v can be taken from 
Levenfish and Yakovlev [42] : 




(A.l) 



V — 



A(")(T) 



1.456 



0.157 



+ 



1.764 



T 



(A.2) 
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In the case of triplet-state neutron pairing the effective gap should be intro- 
duced which is defined by Eq. (56). The function v is now given by (see Ref. 
[30]): 

• - - VT^r fo.7893 + H!^), r ^ ^ . (A.3) 



Gnedin and Yakovlev [26] calculated the function $«(!') in a wide range of 
V (for the problem of thermal conductivity). These authors fitted ^ai^) by a 
simple analytical formula which is correct at any v and satisfies the asymptotes 
$a = \/27n; e"'' at v ^ +oo: 

1/2 



a 



0.9443 + a/(0.0557)2 + (0.1886^;)2 



exp ( 1.753- ^(1.753)2 + 172 

(A.4) 



Calculation and fit errors do not exceed 2.6%. 

B Fits to Fi and F{ 

B.l Symmetric nuclear matter 



In the symmetric nuclear matter S — 0. The plots of functions Fi and F[ 
versus the wave number fepsym — (37r2?T,b/2)^/^ are given by Jackson et al. [33] 
for the model interactions BJ v6 and Rcid v6 (see Section 4.2 for details). 
We have approximated these functions by simple analytical formulas in the 
interval 1.2 fm~^ < kpsym < 2.0 fm~^. For the BJ v6 interaction we obtain 

Fi(nb, 0) = -0.6854 + 0.6724 fcpsym - 0.5180 {kp^y^Y , (B.l) 
F;(nb, 0) = 1.723 - 1.520 A;Fsym + 0.03498 {kp^yj^ , (B.2) 

while for the Reid v6 interaction 

Fi(nb, 0) = 1.034 - IMQkFsym + 0.5455 {k^^yj^ , (B.3) 
F[{nb, 0) = 0.6973 + 0.1403 /cpsym - 0.5303 {kpsymf ■ (B.4) 

B.2 Pure neutron matter 



In the pure neutron matter 6=1. The plots of the function Fi versus the wave 
number of the pure neutron matter /cppure = (37r2nb)^^'^ are given by Jackson 
et al. [33] (see also Section 4.2). 
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The fit for tlie BJ v6 interaction lias tlie form 

FiK, 1) = 0.1473 + 0.7372 h^,,,, - 1.0414 (h^^,,)^ + 0.1958 {h^^^.f , (B.5) 
while for the Reid v6 we obtain 

Fi(nb, 1) = -0.2729+1.5545 fepp^re- 1-3225 {kFp^,,f +0.2393 (fcppure)^ • (B.6) 

The fitting formulas (B.5) and (B.6) correctly describe the function -Fi(r2b, 1) 
in the interval 0.75 fm~ < fcppure ^ 3.0 fm~ . 
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